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fi 要 : [目的 /意义 ] 为 实现 玉米 的 干旱 胁迫 等 作物 生长 状态 的 无 损 监 测 ， 本 研究 探索 一 种 基于 高 光谱 技术 的 干 
旱 胁迫 下 玉米 叶片 叶绿素 含量 与 含水 量 无 损 检 测 方法 。[ 方 法 ] 首先 使 用 高 光谱 相机 采集 不 同 干旱 胁迫 程度 的 苗 划 
玉米 叶片 图 像 ， 并 使 用 图 像 处 理 技术 提取 叶肉 部 分 平均 光谱 。 通 过 系统 性 地 分 析 不 同 特征 波长 提取 方法 、 机 器 学 
习 回 归 模 型 对 叶绿素 含量 和 含水 量 预测 性 能 的 影响 ， 分 别 建立 最 优 叶 绿 素 含量 和 含水 量 反 演 模 型 ， 并 探究 构建 可 
用 于 叶绿素 含量 和 含水 量 反 演 的 植被 系数 并 评估 其 反 演 能 力 。[ 结 果 和 讨论 ] 结合 逐步 回归 (Stepwise Regression, 
SR) 特征 提取 与 Stacking 回归 可 获得 最 优 叶绿素 含量 预测 效果 ( 尼 为 0.878， 均 方 根 误差 为 0.317 mg/g); 结合 连续 投 
影 算法 (Successive Projections Algorithm, SPA) 特征 提取 与 Stacking 回归 可 获得 最 优 含水 量 预测 效果 CRO 
0.859, RMSE 为 3.75%); 新 构建 的 归 一 化 差分 植被 指数 LC Ry Reso) / (Ry ot Rss) | 和 比值 系数 (RyRy) 分 别 
对 叶绿素 含量 和 含水 量 反 演 精度 最 高 且 显 著 高 于 传统 植被 系数 ， 尼 分 别 为 0.803 和 0.827， 均 方 根 误差 分 别 为 
0.403 mg/g 和 3.28%。 [结论 ] 本 研究 构建 的 基于 高 光谱 信息 的 反 演 模型 与 植被 系数 可 以 实现 玉米 叶片 叶绿素 含量 
与 含水 量 的 精确 、 无 损 检 测 ， 可 为 玉米 生长 状态 实时 监测 提供 理论 依据 和 技术 支持 。 
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1 引言 要 生理 指标 ， 可 以 反映 干旱 胁迫 程度 。 因 此 ， 对 玉 
米 叶片 叶绿素 含量 与 含水 量 的 快速 、 无 损 检 测 有 户 
向 作 惫 品质 与 产量 中， 快速 诊 肠 干 里 胁 和 程度 是 及 “实现 靖 准 农业 中 玉米 农田 大 范围 的 干旱 胁迫 在 线 监 


时 采取 调控 措施 以 缓解 干旱 胁迫 的 前 提 。 虽 然 传 统 Vo EARE a a 
的 实地 调研 作物 外 观 形态 以 及 检测 生物 学 指标 等 方 AGE 


玉米 (Zea mays L.) 苗 期 的 干旱 胁迫 会 严重 影 


法 可 以 准确 评估 作物 生长 状态 ， 但 需要 大 量 人 力 物 。“” 详 技 术 和 成 像 技术 结合 ， 使 其 具有 广 波段 荐 商 、 全 
力 ， 不 适用 于 大 范围 的 干旱 胁迫 在 线 监测 。 玉 ”空间 信息 和 图 谱 合 一 的 技术 优势 ， 可 用 于 快速 、 无 
米 叶片 叶绿素 含量 与 仿 水 量 是 玉米 生长 过 程 中 的 重 。 损 检测 物体 的 外 部 特征 和 内 部 成 分 。 目 前 ， 国 内 外 
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开展 了 大 量 基于 高 光谱 成 像 技术 的 作物 生理 信息 反 
演 研 究 。 部 分 研究 在 实验 室 中 采集 叶片 的 高 光谱 图 
像 ， 通 过 提取 叶片 部 分 的 平均 光谱 ， 经 预 处 理 并 提 
取 特 征 后 建立 生理 信息 反 演 模型 。 例 如 ， 孙 俊 等 ” 
比较 了 不 同 预 处 理 方法 ， 并 用 偏 最 小 二 乘 回 归 
(Partial Least Squares Regression，PLSR) 建立 了 生 
菜 叶 片 中 氮 素 含量 的 检测 模型 ; Gao 等 “使 用 相关 
性 分 析 与 随机 蛙 跳 算法 提取 特征 波段 ， 随 后 用 
PLSR 建立 玉米 叶绿素 含量 估 测 模型 。 部 分 研究 使 
用 便携 式 、 和 车载、 无 人 机 等 高 光谱 平台 ,采集 植 
株 、 大 田 尺 度 的 高 光谱 图 像 ， 用 于 建立 基于 高 光谱 
图 像 的 群体 营养 成 分 预测 模型 。 例 如 ，Tan 等 利 
用 无 人 机 平台 拍摄 了 毛 胁 迫 下 小 麦 不 同 生 长 时 期 的 
图 像 ， 通 过 融合 文 持 向 量 回归 与 虫 群 优 化 算法 提取 
光谱 特征 ， 并 基于 PLSR 构建 了 氮 含 量 预 测 模型 ; 
Shu 等 ”构建 并 比较 了 多 种 植被 系数 用 于 判断 玉米 
植株 的 水 分 状况 。 但 现 有 方法 大 多 采用 特定 预 处 理 
和 建 模 方法 构建 作物 生理 信息 反 演 模型 ， 而 对 不 同 
预 处 理 、 特 征 提取 、 建 模 方法 及 其 组 合 在 生理 信息 
反 演 中 的 性 能 进行 系统 、 全 面 的 对 比分 析 研 究 较 
少 ， 尤 其 对 于 在 干旱 胁迫 下 检测 玉米 苗 期 生理 信息 
反 演 方法 的 系统 性 分 析 研 究 则 更 少 。 

本 研究 以 干旱 胁迫 下 的 玉米 植株 为 研究 对 象 ， 
通过 系统 性 地 分 析 不 同 预 处 理 方 法 、 特 征 波段 提取 
方法 、 回 归 模 型 对 含水 量 、 叶 绿 素 含量 预测 性 能 的 
影响 ， 建 立 性 能 最 优 的 基于 高 光 成 像 技 术 的 叶片 叶 
绿 素 含 量 和 含水 量 反 演 模 型 ; 构建 可 用 于 叶绿素 和 
含水 量 反 演 的 植被 系数 ， 并 评估 其 反 演 能 力 。 本 人 研 
究 可 为 玉米 生长 信息 的 快速 获取 提供 技术 支持 ， 以 
期 实现 对 玉米 种 植 的 精细 化 管理 ， 对 提高 玉米 产量 
和 品质 具有 重要 意义 。 


2 材料 与 方法 


2.1 干旱 胁迫 玉米 叶片 样本 制备 


实验 于 2022 年 9 月 在 广西 大 学 农学 院 温 室 大 棚 
进行 。 玉 米 品 种 为 “早生 糯 玉 808”。 在 玉米 苗 期 
(4 片 叶 ) 对 玉米 植株 进行 干旱 胁迫 处 理 。 干 旱 胁迫 
实验 设置 正常 水 处 理 [田间 持 水 量 的 (75+2) %, 
CK], EFF | 田间 持 水 量 的 (6042) %, WIJ, 


中 度 干 旱 [田间 持 水 量 的 (45+2) %，W2]、 重 度 
干旱 [田间 持 水 量 的 (30+2) %, W3] 共 4 个 干旱 
梯度 ,每 组 12 株 ,共计 48 株 。 干旱 胁迫 全 部 在 大 
棚 内 进行 ， 以 防止 降雨 对 土壤 含水 量 的 影响 。 在 胁 
迫 第 3、6、9 天 的 每 天 9 : 00 am 对 玉米 植株 进行 共 
3 次 采样 操作 。 每 次 采样 从 4 组 处 理 中 各 选取 4 株 玉 
米 植株 ， 采 集 每 盆 玉 米 顶 端的 3 片 叶片 ， 用 剪刀 将 
每 片 叶片 剪 成 2 片 ， 随 后 立即 放 和 人 标本 袋 中 带 回 实 
验 室 。 共 采集 玉米 叶片 样本 总 数 3X4X4X3X2= 
288 个 。 


2.2 高 光谱 图 像 采集 


本 研究 使 用 的 高 光谱 成 像 系 统 由 光源 平行 放置 
的 两 个 高 光谱 相机 Specim FX10 和 Specim FX17 
(Spectral Imaging Ltd., Oulu, Finland), 、 相 机 推 扫 
运动 控制 系统 和 计算 机 组 成 (图 1)。 高 光谱 相机 的 
主要 技术 参数 见 表 1。 在 采集 高 光谱 图 像 前 ， 使 用 
标准 反射 率 的 白板 与 黑板 进行 测量 标定 。 图 像 采 和 集 
在 昏暗 环境 下 (环境 照度 <5 lux) 进行 ， 并 在 采集 
过 程 中 保持 环境 光 条 件 一 致 。 通 过 将 每 个 叶片 样 
本 置 于 黑色 背景 的 样本 架 中 ， 控 制 电机 人 匀速 运动 
带动 相机 推 扫 以 采集 高 光谱 图 人像。 分 别提 取 
FX10 与 FX17 相机 采集 的 400~1000 nm 波段 与 
1000~1700 nm 波段 ,经 过 后 处 理 配 准 倒 加 为 一 张 
高 光谱 图 像 ， 共 生成 玉米 叶片 光谱 照片 288 张 。 

表 1 高 光谱 相机 主要 技术 参数 


Table 1 Specs of hyperspectral cameras 


ee 相机 型 号 
技术 参数 
FX10 FX17 
光谱 范围 mm 400~ 1000 900~1700 
光谱 分 辨 率 /nm 5.5 8 
狭 颖 宽度 /um 30 30 
空间 像素 数 1024 640 


像素 大 小 /um 8X8 


2.3 玉米 叶片 叶绿素 含量 与 含水 量 测定 


高 光谱 图 像 采 集 完毕 后 ,使 用 标准 方法 测定 叶 
片 中 叶绿素 a 和 叶绿素 b 含 量 ”。 将 每 片 叶片 避 开 
叶脉 部 分 迅速 剪 取 0.1 g， 揭 碎 后 放 和 人 含有 25 mL 的 
丙酮 和 无 水 乙醇 (体积 比 2 : 1) 茶 取 液 的 试剂 管 
中 ， 置 于 阴暗 柜子 中 24h， 使 叶绿素 萃取 分 离 后 取 
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“ey, mg/g。 计 算 方 法 如 公式 (4) 所 示 。 
— _VxC 

~ Ton i 
<0 必 区 | S Wee 其 中 ，7 为 某 取 液体 积 ，L; m 为 叶片 质量 ,mg。 
人 电机 控制 器 yE 
HEH WY 将 叶片 去 除 叶 脉 ， 使 用 电子 天 平 称 取 叶 片 鲜 重 

wea ~ = ` ETN y 
系统 质量 (W)， 再 将 叶片 放 入 落 馏 水 中 2h， 使 其 完全 


上 


(a) 光 谱 相 机 与 光源  (b) 推 扫 控 制 器 与 样本 放置 处 
图 1 用 于 采集 玉米 叶片 光谱 的 高 光谱 成 像 系统 
Fig. 1 The hyperspectral imaging System used to collect maize 


leaf spectrum 


出 。 待 叶片 完全 发 白 后 ， 将 萃取 液 转 移 至 比 色 严 
中 ， 放 入 分 光 光 度 计 (UV-1800， 
ration，Kyoto，Japan) ， 检 测 萃取 液 在 663 和 645 nm 
处 的 吸光 度 。 通 过 公式 (1) ~ 公式 (3) 计算 叶 绿 


素 含量 : 


Shimadzu Corpo- 


C, = 12.724 6 — 2.59A eas (1) 
C, = 22.884 s — 4.67A os (2) 
C=C, + C, =20.294 a; + 8.054 ges (3) 


其 中 ，4 FA gas IIJ ARTE 663 和 645 nm 
处 的 吸光 度 ; C, 和 C, 分 别 为 溶液 中 叶绿素 a 和 叶 绿 
素 b 的 含量 ，mg/L; C 为 溶液 的 叶绿素 含量 ，mg/L。 
溶液 中 叶绿素 含量 需要 转换 为 叶片 叶绿素 含量 Chl, 


膨胀 ， 取 出 叶片 ， 将 叶片 擦 干 ， 再 次 称 取 叶片 饱和 
质量 (WW,)， 最 后 将 叶片 放 入 干燥 箱 中 ，105°C 杀 青 
0.5 h，80°C 便 温 烘 烤 至 叶片 恒 重 。 此 时 的 质量 为 叶 
片 干 重 质量 (Wp) IMF (相对 ) 含水 量 (Liquid 
Water Content, LWC) 计算 方法 如 公式 (5) 所 示 。 


W.-W, 
ma (5) 


LWC = 


2.4 数据 处 理 方法 


采集 的 高 光谱 图 像 中 叶肉 、 叶 脉 和 背景 部 分 的 
光谱 存在 明显 的 差异 (图 2)。 其 中 背景 部 分 光谱 整 
体 反 射 率 低 于 0.1。 叶 肉 部 分 的 光谱 曲线 呈现 典型 
的 植物 特性 ， 具 有 绿 峰 、 红 和 谷 、 红 边 等 特征 。 叶 脉 
部 分 在 400~-1300 nm 的 光谱 反射 率 整 体高 于 叶肉 光 
谱 反 射 率 ; 而 在 1300~1700 nm 的 光谱 反射 率 低 于 
叶肉 。 由 于 叶片 叶脉 部 分 叶绿素 含量 很 低 ， 因 此 本 
研究 对 高 光谱 图 像 进 行 分 害处 理 以 去 除 叶 脉 篆 
景 ， 只 保留 叶肉 部 分 用 于 提取 光谱 信息 。 


(a) 原 始 叶片 图 像 与 像素 类 型 


由 于 背景 、 叶 肉 和 叶脉 的 光谱 反射 率 差 异 主要 
集中 在 550 和 1100 nm 波段 ， 因 此 在 提取 以 上 两 波 
段 的 反射 率 图 像 时 使 用 大 津 法 分 制 叶脉 PA AF 
景 , 保留 叶肉 区 域 图 像 (图 2 (c) )， 并 提取 区 域 
内 像素 ， 同 时 计算 叶肉 平均 光谱 。 

本 研究 采用 “ 预 处 理 - 特 征 提 取 - 机 器 学 习 反 演 ” 
的 算法 流程 框架 (图 3)。 首 先 ， 对 光谱 进行 预 处 理 


0-200 600 800 1000 1200 1400 1600 
波长 /nm 

(b) 不 同 像素 类 型 对 应 的 光谱 
图 2 叶片 各 区 域 光 谱 差 异 以 及 叶肉 部 分 分 割 结 


Fig. 2 The spectral differences for each leaf segment and the segmentation results for leaf lamina area 


(c) 叶 肉 部 分 分 割 效 果 图 


以 消除 环境 、 仪 器 等 因素 造成 的 噪声 ， 突 出 光谱 中 
有 用 信息 。 本 研究 分 析 的 预 处 理 方 法 包括 卷 积 平滑 
法 ” (Savizky-Golay, SG). 4 JG # Ht BIE” 
(Multiplicative Scatter Correction, MSC). 28 8th 
准 化 "(Standard Normal Variate，SNV) 、 一 阶 导 
数 (First Derivative, FD) 和 二 阶 导 数 (Second 
Derivative, SD). 
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其 次 ， 对 预 处 理 后 的 光谱 进行 特征 提取 与 数据 
降 维 以 减少 光谱 曲线 中 的 元 余 信 息 ， 防 止 过 拟 合 。 
光谱 数据 降 维 与 特征 提取 方法 包括 特征 波长 提取 和 
植被 系数 法 ， 特 征 波长 提取 包括 连续 投影 算法 '” 
(Successive Projections Algorithm, SPA) 、Pearson 
相关 系数 法 、 随 机 森林 算法 (Random Forest, RF) 
和 逐步 回归 分 析 (Stepwise Regression，SR )。 植 被 
系数 (Vegetation Index, VI) 是 某 些 特定 波段 反射 
率 的 线性 或 非 线性 组 合 。 本 研究 分 析 传 统 植被 系数 
( 表 2) 对 叶绿素 含量 和 含水 量 的 反 演 能 力 ， 并 对 几 
种 党 用 植被 系数 〈 表 3) 的 波段 进行 重新 调整 以 构 
建新 的 植被 系数 ， 包 括 差 值 系数 (Difference In- 
dex，DI) 、 比 值 系数 (Ratio Index, RI) 和 归 一 化 
差 值 植被 系数 (Normalized Difference Vegetation In- 
dex，NDVI)， 并 评估 其 反 演 能 

表 2 常用 传统 植被 系数 


Table 2 Traditional vegetation indices 


、 参考 

植被 名 称 及 计算 公式 描述 文献 

GIFR 554/Roo7 (6) [13] 
CL_730 = Rgoo/R730—1.0 (7) [14] 
CI 709 = Rego /Ry99—1-0 (8) 叶绿素 [4] 
Chl_ green = Raoo/Rss0—1.0 (9) ”会 量 相关 的 ”[15] 


NDRE =(Ro-R0)/ (Roo + Roy) (10) ”植被 系数 。 [16] 


RGVI = R;so/R67 (11) [17] 
MTCI=(Ry54 — Rao )/ (Rr + Rei) (12) [18] 
NDWI = (Rgs R1240 )/ (Rgoo + Rizo) (13) [19] 

MST =R 600R s20 (14) [20] 

i 含水 量 y 

ANDVI = ( Rooy—Rego )/ (Roo + Rego) (15) 、 n [21] 
相关 的 植被 

WI = Ro00/Ro70 (16) 系数 es 

SRWI = Rg20/R 1200 (17) [23 | 
NDI = ( Rg20=R 1649 )/ (Rao + Risso) (18) [21] 


TE: R; 为 波长 inm 对 应 的 光谱 反射 率 


RI 植被 系数 名 称 及 公式 


Table 3 Vegetation index names and formulas 


植被 系数 名 称 及 计算 公式 参考 文献 
DI =R,/R, (19) [24] 
RI=R,/R, (20) [24 | 
NDVI =(R; + R/R, + R,) (21) [24] 


TE: R,、R 分别 代表 波长 i 和 j nm 对 应 的 反射 率 


随后 ， 利 用 机 器 学 习 技术 构建 叶绿素 含量 和 合 
水 量 反 演 模型 。 机 器 学 习 方 法 包括 PLSR、 神 经 网 
络 (Artificial Neural Network，ANN) 、 支 持 向 量 回 
J (Support Vector Regression，SVR) 、 玉 近邻 (K- 
nearest Neighbor，KNN) 和 Stacking。 

最 后 ， 通 过 系统 比较 不 同 预 处 理 方法 、 数 据 降 
维 方法 与 机 器 学 习 模 型 的 性 能 ， 选 择 性 能 最 优 的 叶 
绿 素 含 量 和 含水 量 反 演 模型 。 本 研究 随机 将 70% 的 
光谱 样本 数据 (202 个 样本 ) 作为 训练 数据 集 ， 用 
于 比较 不 同 预 处 理 和 特征 提取 方法 ， 并 训练 反 演 模 
型 ;将 剩余 30% 的 光谱 样本 数据 (86 个 样本 ) 作为 
测试 数据 集 ， 用 于 评估 反 演 模型 性 能 。 使 用 的 模型 
性 能 评价 指标 包括 决定 系数 R 和 均 方 根 误差 
(Root Mean Squared Error, RMSE). 


两 光谱 图 像 


分 割 


亨 优 似 处 理 后 的 高 光谱 数据 


| 
ISEA GI / [ aame ooo 广 一 一 


¥. Y. x 
[nasa ] 相关 系数 法 | | menem | | 随机 森林 算法 | | esun 
k 4 
| mwenn j 
+ ¥ t } 
线性 回归 ANN SVR KNN Stacking 
I 


特征 波段 与 机 器 学 习 结合 
的 模型 (16 种 ) 


评价 指标 R* 与 RMSE 
图 3 叶绿素 含量 和 含水 量 反 演 研究 整体 流程 
Fig. 3 The overall process for the study of the estimation of 


leaf chlorophyll content and water content 


3 ”结果 与 讨论 


3.1 不 同 程度 干旱 胁迫 下 的 玉米 叶片 光谱 特 
征 与 变化 趋势 分 析 


受 胁迫 与 未 受 胁迫 叶片 表现 出 不 同 的 光谱 特征 
(图 4) ， 即 随 着 胁迫 的 加 重 ， 叶 片 反 射 率 在 波长 
400~1700 nm 范围 内 随 干旱 胁迫 程度 逐步 递增 ， 但 
光谱 曲线 总 体 趋势 一 致 。 
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该 现象 的 植物 生理 学 原因 可 分 波段 解释 : 
400~750 nm 波段 反射 率 受 叶绿素 吸收 作用 影响 ， 
干旱 胁迫 会 导致 玉米 叶片 叶绿素 分 解 ， 使 叶片 的 叶 
绿 素 含 量 降低 ; 750~1350 nm 波段 的 反射 率 主要 受 
叶片 内 部 结构 的 影响 , 干旱 胁迫 会 使 叶片 内 部 结构 
发 生 损伤 ; 1350~1700 nm 波段 的 反射 率 受 水 的 吸 
收 作用 影响 ,干旱 胁迫 会 降低 叶片 的 含水 量 。 


0.0 
400 600 800 1000 1200 1400 1600 1800 
波长 /nm 


图 4 不 同 干旱 胁迫 下 玉米 叶片 光谱 曲线 


Fig. 4 Spectral profiles of maize leaves under different 


drought stresses 


3.2 不 同 光谱 预 处 理 方法 的 性 能 比较 


将 SG、MSC、SNV、FD、SD、SG+MSC、 
SG+SNV,, SG+ED 和 SG+SD 共 9 种 算法 用 于 光谱 数 
据 预 处 理 (图 5)， 可见 SG 平滑 后 的 光谱 曲线 与 原 
始 光 谱 差 异 不 明显 。 这 是 由 于 提取 光谱 为 叶肉 部 分 
的 平均 光谱 ， 噪声 较 少 。 而 经 过 MSC 和 SNYV 归 一 
化 处 理 后 的 光谱 分 布 变 得 集中 。 经 过 FD 和 SD 差分 
处 理 后 的 光谱 吸收 峰 的 数量 明显 增加 。 

为 评估 不 同 预 处 理 方法 的 性 能 ， 将 原始 光谱 数 
据 和 其 他 9 种 预 处 理 后 的 光谱 数据 作为 PLSR 的 输 
入 ,叶绿素 含 量 和 含水 量 作为 PLSR 的 输出 ， 通 过 
5 折 交 叉 验 证 确定 最 佳 主 成 分 数 参 数 ， 并 使 用 验证 
集 评估 预测 性 能 ， 进 而 确定 最 优 的 光谱 预 处 理 方法 
( 表 4)。SG+MSC 预 处 理 对 叶绿素 含量 预测 性 能 最 
FE; 原始 光谱 数据 对 含水 量 预测 性 能 最 佳 。 因 此 ， 
在 后 续 反 演 模型 的 建立 中 ， 采 用 SG+MSC 预 处 理光 
谱 数 据 用 于 叶绿素 反 演 ;采用 原始 光谱 数据 用 于 含 
水 量 反 演 。 


at 0.3 ay 
03 E p gI 
me ri Eos i 
0. -of 25 
è -01 
00 “15 10 
s00 1200 1600 400 s00 1200 1600 -2900 800 1200 1600 400 30) 1200 1600 -0200 S00 1200 1600 
波长 /nm 波长 /nm 波长 /nm 波长 /nm 波长 nm 
(a) 原 始 光 谱 (b)MSC 处 理 (c)SNV 处 理 (d)FD 处 理 (e)SD 处 理 
10 x10 
10 
_ 0S 
WY R00 
1 
-0.5 
-10 
BOO 1200 1600 -2900 200 1200 1600 400 300 1200 1600 400 300 1200 160) 
波长 /nm 波长 :nm 波长 /nm 波长 mm 波长 mnm 
(PASGA (g)SG+MSC 处 理 (h)SG+SNV 处 理 (i)SG+FD 处 理 (j)SG+SD 处 理 


图 $ 不 同 预 处 理 后 的 玉米 叶片 反射 光谱 曲线 图 


Fig. 5 Plots of reflectance spectra of corn leaves using different preprocessing algorithms 


3.3 基于 特征 波长 的 玉米 叶片 叶绿素 含量 和 
含水 量 反 演 模型 构建 
3.3.1 特征 波长 提取 结果 分 析 
不 同 特征 波长 提取 方法 第 选 获 得 的 波段 如 图 6 
所 示 。 其 中 SPA 算 法 以 RMSE 最 小 值 对 应 变量 数 确 


定 提取 的 特征 波长 数量 。 叶 绿 素 含量 反 演 提取 的 15 
个 特征 波长 分 布 于 整个 波段 范围 内 (图 6 (a) ); 
含水 量 反 演 提取 的 16 个 特征 波长 主要 集中 在 大 于 
900 nm 的 波长 范围 (图 7 (a) )。Pearson 相关 系数 
法 计算 各 波长 反射 率 与 叶绿素 含量 和 含水 量 的 相关 
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表 4 不 同 预 处 理 后 光谱 与 原始 光谱 对 叶绿素 含量 和 含水 量 
预测 性 能 
Table 4 Inversion performance of different pre-treated and orig- 


inal Spectra for chlorophyll content and water content 


预 处 理 ER 训练 集 验证 集 
反 演 目标 ` p 
方法 分 数 R? RMSE, RÌ RMSE, 
Te 10 0.857 0.307 0.839 0.312 
SG 10 0.866 0.298 0.832 0.318 
MSC 11 0.867 0.297 0.831 0.312 
SNV 12 0.873 0.289 0.823 0.326 
叶绿素 含 FD 6 0.908 0.247 0.838 0.312 
Ht/(mg-g |) SD 4 0.908 0.246 0.809 0.339 
SG+MSC 11 0.857 0.307 0.839 0.311 
SG+SNV 12 0.872 0.291 0.825 0.324 
SG+FD 6 0.884 0.277 0.830 0.320 
SG+SD 4 0.898 0.260 0.807 0.341 
无 11 0.877 3.05 0.843 3.57 
SG 11 0.877 3.09 0.843 3.58 
MSC 13 0.856 3.34 0.796 4.07 
SNV 14 0.863 3.26 0.702 4.92 
FD 5 0.849 3.42 0.624 5.53 
含水 量 /% 
SD 4 0.868 3.20 0.710 4.85 
SG+MSC 13 0.855 3.35 0.796 4.07 
SG+SNV 14 0.861 3.28 0.706 4.89 
SG+FD 5 0.839 3.53 0.632 5.47 
SG+SD 4 0.862 3.27 0.734 4.65 
一 nae 0.45 一 jaa 
© 特征 波长 点 0.40 ”特征 波长 点 
0.35 
St 0.30 
aR 0.25 
1X 0.20 
O15 
0.10 
0.05 R r í 
400 600 800 1000 1200 1400 1600 400 600 800 1000 1200 1400 1600 
波长 /nm 波长 /nm 
(a)SPA 算法 (b) 相 关系 数 法 


系数 ， 根 据 其 绝对 值 从 大 到 小 对 波长 排序 。 随 后 依 
次 选取 不 同 数 量 波长 与 线性 回归 方法 结合 ， 分 别 建 
立 叶绿素 含量 和 含水 量 的 反 演 模型 以 评 佑 最 优 特征 
波长 数量 。 最 终 选取 相关 系数 绝对 值 通过 将 相关 系 
数 绝对 值 从 大 到 小 排序 ， 选 取 前 25 个 波长 作为 叶 
绿 素 含量 反 演 的 特征 波长 (图 6 b) ); 选取 相关 系 
数 绝对 值 前 35 的 波长 作为 含水 量 反 演 的 特征 波长 
(图 7 (b) )。RF 算 法 计算 每 个 波长 对 叶绿素 含量 、 
含水 量 反 演 的 ， 表征 变量 重要 性 的 基尼 系数 。 
对 于 叶绿素 含量 以 及 含水 量 反 演 ， 基 尼 系 数 较 高 
的 波长 分 别 集中 在 400~800 nm 波段 以 及 800~ 
1350 nm 波段。 根据 基尼 系数 从 大 到 小 对 波长 排 
序 ， 随 后 依次 选取 不 同 数量 波长 与 线性 回归 方法 结 
合 , 分 别 建立 叶片 叶绿素 含量 和 含水 量 的 反 演 模型 
以 评估 最 优 特征 波长 数量 。 最 终 选 取 基 尼 系 数 从 大 
到 小 前 20 的 波长 作为 叶绿素 含量 反 演 的 特征 波长 
(图 6 C) ); 从 大 到 小 选取 前 25 个 基尼 系数 对 应 
的 波长 作为 含水 量 反 演 的 特征 波长 (图 7 (c) )。 
SR 算法 为 叶绿素 含量 反 演 选择 了 9 个 特征 波长 ， 分 
布 于 整个 波段 范围 内 (图 6 d) ); 为 含水 量 的 反 
演 选 择 了 13 个 特征 波长 ， 集 中 在 1300 nm 波长 附近 
(图 7 (d) )。 


一 一 kah 
ee。 特征 波长 点 


400 600 800 1000 1200 1400 1600 


400 600 800 1000 1200 1400 1600 


波长 mm 波长 /nm 
(c)RF 算 法 (d)SR 算 法 


图 6 各 算法 提取 用 于 叶绿素 含量 反 演 的 特征 波长 分 布 


Fig. 6 Distributions of the characteristic wavebands extracted for chlorophyll content inversion by algorithms 


3.3.2 ”基于 特征 波长 的 反 演 模型 性 能 分 析 

为 比较 不 同 特征 波段 提取 方法 和 回归 方法 对 叶 
绿 素 含 量 和 含水 量 反 演 性 能 的 影响 ， 将 基于 SPA、 
相关 系数 、RF、SR 这 4 种 方法 提取 的 特征 波长 以 
及 全 波长 共 5 组 变量 , 分别 作 为 ANN、SVR、KNN 
和 Stacking 共 4 种 回归 方法 的 输入 。 分 别 以 叶绿素 


含量 和 含水 量 作为 输出 ， 各 建立 20 种 模型 。 其 中 
Stacking 的 初始 学 习 央 采用 ANN 、SVR、KNN 和 线 
性 回归 器 进行 预测 。 通 过 $ 折 交叉 验证 确定 最 佳 参 
数 ， 使 用 测试 集 评估 模型 预测 性 能 ， 进 而 分 别 确定 
性 能 最 优 的 叶绿素 含量 和 含水 量 的 反 演 模型 方法 。 

对 于 叶绿素 含量 的 预测 ， 不 同 模型 性 能 如 表 5 
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图 7 各 算法 提取 用 于 含水 量 反 演 的 特征 波长 分 布 


Fig.7 Distributions of the characteristic wavebands extracted for water content inversion by each algorithm 


表 5 不 同 模 型 对 玉米 叶片 叶绿素 含量 的 预测 性 能 


Table 5 Inversion results of different models for chlorophyll 


content of maize leaves 


T 训练 集 测试 集 
特征 波长 ”回归 
提取 方法 方法 Bere ga IE 
(mg:g') ¥ (mg: g") 
ANN 0.849 0.292 0.812 0.394 
SVR 0.854 0.288 0.865 0.334 
a KNN 0.831 0.309 0.766 0.439 
Stacking 0.870 0.272 0.849 0.353 
ANN 0.898 0.241 | 0.816 加 0.390 
SVR 0.845 0.297 0.806 0.400 
SPA 
KNN 0.844 0.298 0.750 0.455 
Stacking 0.869 0.273 0.840 0.364 
ANN | 0.833 | 0.308 o 0.773 0.433. | 
SVR 0.752 0.375 0.741 0.462 
相关 系数 
KNN 0.596 0.479 0.587 0.584 
Stacking 0.761 0.368 0.775 0.431 
SVR 0.764 0.366 0.735 0.468 
KNN 0.615 0.467 0.594 0.579 
Stacking 0.791 0.343 0.778 0.428 
SVR 0.828 0.312 0.849 0.353 
= KNN 0.810 0.328 0.796 0.410 
Stacking 0.855 0.287 0.878 0.317 


所 示 。 通 过 比较 


可 知 ， 结 合 SR 特征 提取 和 Stack- 


ing 回归 的 模型 为 反 演 叶绿素 含量 性 能 最 优 的 模 
型 ， 测 试 集 决定 系数 及 为 0.878， 测 试 集 RMSE ,为 
0.317 mg/g (图 8) 与 使 用 全 波段 Stacking 模型 相 


lt, 454 SR FA 


FE 提 取 和 Stacking 回归 的 模型 不 仅 


在 测试 集 上 使 决定 系数 她 提高 2.90%，RMSE, 降 低 


0.036 mg/g， 同 时 使 模型 的 输入 变量 数 降 为 9， 减 少 
99.31%。 由 此 可 以 证 明 ，SR 特征 提取 算法 能 够 在 
叶绿素 含量 预测 中 有 效 地 提取 特征 波长 ERIR 
信息 ， 从 而 提高 预测 性 能 。 

对 于 含水 量 的 预测 ， 不 同 模型 性 能 如 表 6 所 
示 。 通 过 比较 ， 并 考虑 模型 输入 参数 量 ， 可 以 确定 
以 ANN 和 SVR 为 回归 方法 时 ， 结 合 SPA 特征 提取 
和 Stacking 回归 的 模型 为 反 演 叶绿素 含量 性 能 最 优 
的 模型 ， 测 试 集 决 定 系数 有 尼 为 0.878， 测 试 集 
RMSE, 为 0.317 mg/g (图 9)。 与 使 用 全 波段 Stack- 
ing 模 型 相 比 ， 结 合 SPA 特征 提取 和 Stacking 回归 的 
模型 不 仅 在 测试 集 上 决定 系数 已 提高 0.2%、 
RMSE, 降 低 0.03%， 同 时 使 模型 的 输入 变量 数 降 为 
16, 减少 98.77%。 由 此 可 以 证 明 ，SPA 特征 提取 算 
法 能 够 在 含水 量 预测 中 有 效 地 提取 特征 波长 ， 去 除 
匈 余 信息 ， 从 而 提高 预测 性 能 。 
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图 8 使 用 SR-Stacking 模 型 对 测试 集 数据 的 叶绿素 含量 预测 
结果 散 点 图 


Fig.8 Scatter plot of chlorophyll content prediction results us- 


ing the SR-Stacking model on the testing dataset 
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表 6 不 同 模型 对 玉米 叶片 含水 量 的 预测 性 能 
Table 6 Inversion performance of different models for the wa- 


ter content of maize leaves 


寺 征 波长 训练 集 测试 集 
it 可 归 方 法 - - 
提取 方法 R? RMSE/% Ri  RMSE,/% 
ANN 0.767 3.99 0.835 4.06 
SVR 0.709 447 0.799 4.47 
无 
KNN 0.777 3.92 0.804 4.41 
Stacking 0.792 3.78 0.857 3.78 
ANN 0.746 4.19 0.809 4.37 
SVR 0.710 4.47 0.784 4.64 
SPA 
KNN 0.812 3.60 0.827 4.16 
Stacking 0.815 3.58 0.859 3.75 
ANN 0.597 5.28 0.711 5.37 
SVR 0.508 5.83 0.595 6.36 
相关 系数 
KNN 0.643 4.96 0.699 5.48 
Stacking 0.639 4.99 0.724 5.24 
ANN 0.604 5.23 0.699 5.47 
SVR 0.536 5.66 0.539 6.77 
RF 
KNN 0.639 4.99 0.695 5.51 
Stacking 0.633 5.04 0.762 4.87 
ANN 0.753 4.13 0.821 4.22 
SVR 0.719 4.40 — 0.805 4.41 
SR 
KNN 0.764 4.03 0.781 4.67 


Stacking 0.785 


— =x 
R; = 0.8589 
RMSE, = 3.75% 


预测 值 /% 


50 60 70 80 90 100 
真 值 /% 
图 9 使 用 SPA-Stacking 模型 对 测试 集 数据 的 含水 量 预测 结 
果 散 点 图 


Fig. 9 Scatter plot of water content prediction results using the 


SPA-Stacking model on the testing dataset 


3.4 基于 植被 系数 的 玉米 叶片 含水 量 和 时 绿 
素 含 量 反 演 模型 构建 


使 用 预 处 理 后 的 光谱 数据 构建 如 表 2 所 示 的 传 
统 植被 系数 ， 通 过 线性 回归 建立 叶绿素 含量 和 含水 
量 的 反 演 模型 ， 各 传统 植被 系数 的 预测 性 能 如 表 7 
所 示 。 所 有 基于 植被 系数 建立 的 线性 回归 反 演 模型 
尼 都 低 于 0.5， 说 明 其 对 叶绿素 含量 和 含水 量 预测 
性 能 较 差 。 

表 7 不 同 植被 系数 对 叶绿素 含量 和 含水 量 的 预测 性 能 


Table 7 Performance of different vegetation indices in the in- 


version of chlorophyll content and water content 


bats ee 训练 集 测试 集 
反 演 目 标 植被 系数 
Rè RMSE, Rè RMSE, 
GI 0.130 0.704 0.015 0.902 
CI 730 0.281 0.640 0.150 0.838 
CI 709 0.367 0.601 0.292 0.765 
叶绿素 含量 / 
4 Chl green 0.523 0.522 0.438 0.687 
(mg'g ) 
NDRE 0.300 0.632 0.187 0.819 
RGVI 0.126 0.706 0.012 0.903 
MTCI 0.271 0.645 0.208 0.809 
NDWI 0.034 8.25 0.015 10.74 
MSI 0.324 614 0.150 9.51 
hNDVI 0.204 7.01 0.293 8.20 
含水 量 /% 
WI 0.400 5.68 0.338 7.80 
SRWI 0.300 631 0.187 9.13 
NDI 0.126 7.55 0.012 10.77 


因此 ， 探 究 重 新 选择 传统 植被 系数 中 的 波长 ， 
进而 构建 新 的 植被 系数 用 于 叶绿素 含量 和 含水 量 预 
测 的 可 行 性 。 对 于 新 的 植被 系数 构建 ， 首 先 根据 
表 3 中 各 植被 系数 的 计算 公式 ,计算 400~~1700 nm 
范围 内 所 有 这 两 种 波长 组 合 的 差 值 系 数 (Differ- 
ence Index，DI) 、 比 值 系数 (Ratio Index, RI) 和 
归 一 化 植被 系数 (Normalized Difference Vegetation 
Index, NDVI); 其 次 分 析 计 算 获 得 的 DI、RI、 
NDVI 与 叶绿素 含量 和 含水 量 的 相关 性 (图 10); 
对 于 叶绿素 而 言 ， 用 于 计算 植被 系数 的 波段 为 可 见 
光 ; 而 对 于 含水 量 而 言 ， 用 于 计算 植被 指数 的 波段 
为 可 见 光 和 近 红 外 光 。 最 后 选取 相关 系数 绝对 值 最 
大 的 波段 组 合 构建 新 植被 系数 ( 表 8)。 
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(e)NDVI 与 叶绿素 含量 (f)NDVI 与 含水 量 
图 10 新 构建 植被 系数 与 叶绿素 含量 和 含水 量 相关 系数 
热力 图 
Fig. 10 Heat maps of correlation coefficients between the new- 
ly constructed vegetation coefficients and chlorophyll content 
and water content, respectively 
表 8 新 构建 的 植被 系数 对 应 的 相关 系数 最 大 值 及 波长 位 置 
Table 8 The maximum correlation coefficient and wavelength 


position corresponding to the newly constructed vegetation co- 


efficient 
植被 叶绿素 含量 含水 量 
系数 (RI ， ”波长 位 置 (i,j)/mm [RI 波长 位 置 (i,j)/nm 
DI 0.882 (420,558) 0.783 (742,1681) 
RI 0.892 (420,559) 0.898 (400,1171) 
NDVI 0.889 (410,559) 0.890 (410,1348) 


将 新 构建 的 植被 系数 分 别 用 于 建立 叶绿素 含量 
和 含水 量 的 线性 回归 反 演 模型 。 结 果 表 明 ， 当 使 用 
新 构建 的 植被 系数 时 ， 预 测 性 能 相 比 于 使 用 传统 植 
被 系数 的 预测 性 能 具有 明显 提升 。 对 于 叶绿素 含量 
反 演 ， 基 于 归 一 化 植被 系数 NDVI (410, 559) 建 
立 的 回归 模型 预测 精度 最 高 ，R 0.803, RMSE, 
为 0.403 mg/g。 对 于 含水 量 反 演 ， 基 于 比值 植被 系 


数 RI (400, 1171) 建立 的 回归 模型 预测 精度 最 高 ， 
R? J 0.827, RMSE, 43.28% (#9). 


3.5 玉米 叶片 生理 信息 可 视 化 


8.0 mg/g 100% 
6.6 mge 88.0% 
5.2 mg/g 76.0% 
3.8 mg/g 64.0% 
2.4 mg/g 52.0% 
1.0 mg/g 40.0% 


(a) 叶 绿 素 含 量 
All 玉米 叶片 叶绿素 含量 和 含水 量 分 布 热 力图 


Fig. 11 Heatmaps of the distribution of chlorophyll content and 


(b) 含 水 量 


water content of a corn leaf 


本 研究 使 用 SR 提取 的 9 个 特征 波长 和 SR- 
Stacking 模型 对 叶片 叶绿素 含量 进行 可 视 化 研究 ; 
使 用 SPA 提取 的 16 个 特征 波长 和 SPA-Stacking 模 型 
对 叶片 含水 量 进行 可 视 化 研究 。 通 过 对 高 光谱 图 像 
中 叶肉 的 每 个 像素 点 进行 分 析 ， 将 SR 选择 的 9 个 
特征 波长 的 反射 率 作为 SR-Stacking 模 型 的 输入 ， 
输出 每 个 像素 对 应 的 叶绿素 含量 ; 使 用 SPA 选择 的 
16 个 特征 波长 的 反射 率 作 为 SPA-Stacking 模 型 的 输 
入 ， 输 出 每 个 像素 对 应 的 含水 量 ， 并 以 此 构建 叶片 
的 叶绿素 含量 与 含水 量 分 布 图 (图 11)。 结 果 表 明 ， 
对 叶片 生理 参数 的 可 视 化 ， 可 以 更 加 详细 、 直 观 地 
发 现 同 一 叶片 不 同 区 域 之 间 生 理 参数 的 差异 ， 可 用 
于 玉米 的 生长 监测 、 病 害 检 测 、 植 物 育 种 等 研究 
领域 。 


4 结 论 


本 研究 以 玉米 苗 期 叶片 为 对 象 ， 研 究 了 基于 高 
光谱 成 像 搁 术 的 干旱 胁迫 下 叶片 生理 信息 (叶绿素 
含量 和 含水 量 ) 的 快速 无 损 检 测 方法 和 叶片 生理 信 
息 分 布 的 可 视 化 方法 。 首 先 基 于 图 像 处 理 技术 构建 
高 光谱 图 像 分 割 算法 ， 用 于 提取 叶肉 区 域 并 提取 平 
均 反 射 光 谱 ， 排 除 图 像 中 背景 与 叶脉 部 分 的 干扰 ; 
其 次 ， 通 过 系统 地 分 析 不 同 特征 波长 提取 方法 和 机 
顺 学 习 回 归 模 型 对 叶绿素 含量 和 含水 量 的 预测 性 能 
的 影响 ， 分 别 建立 了 最 优 性 能 的 叶绿素 含量 和 含水 


干旱 胁迫 


Vol. 5, No. 3 EMSS : 


自 下 玉米 叶片 叶绿素 含 


含水 量 高 光谱 成 像 反 演 方法 151 


量 的 反 演 模型 ， 最 后 ， 探 究 了 构建 可 用 于 叶绿素 含 


能 力 。 


量 、 仿 水量 的 反 演 的 植被 系数 并 评估 了 其 反 演 
表 9 新 构建 植被 系数 对 叶绿素 含量 和 含水 量 预测 性 能 
Table 9 The inversion performance ofthe newly constructed vegetation coefficients for chlorophyll content and water content 
7 ree 训练 集 测试 集 
反 演 目标 植被 系数 可 归 方 程 
Ri RMSE, R? RMSE, 
DI(420,558) y=28.96x+1.959 (22 0.769 0.363 0.785 0.421 
叶绿素 含量 /mg'g ') RI(420,559) y= -8.691x+10.914 (23 0.774 0.358 0.791 0.415 
NDVI(410, 559) y=15.13x42.268 (24 0.784 0.351 0.803 0.403 
DI( 742, 1681) y=2.014x+1.340 (25 0.731 3.42 0.799 3.53 
含水 量 /% RI(400, 1171) y= —0.9413x+1.5327 (26 0.788 3.04 0.827 3.28 
NDVI(410, 1348) y= —1.526x+0.5752 (27 0.775 3.12 0.807 3.46 
人 研究 结果 表明 ， 将 光谱 技术 用 于 干旱 胁迫 玉米 12]. https://kns.cnki.net/kcms/detail/31.1690.TN.2022103 
Wy Seng sop 1.1514.004.html 
叶片 生理 参数 的 分 析 能 够 实现 生理 参数 的 准确 、 快 ， , 
eee [3] 杨 品 唱 , 卢 俊 玮 . 光谱 技术 在 花生 品质 检测 中 的 应 用 研 
速 检 测 。 对 于 叶绿素 含量 检测 ， 使 用 SG+MSC 预 处 


理光 谱 、SR 提 取 特征 波长 、 建 立 Stacking 模 型 可 获 
得 最 优 性 能 。 对 于 含水 量 检 测 ， 使 用 原始 光谱 ， 
SPA 提取 特征 波长 ， 建 立 Stacking 模 型 可 获得 最 优 
性 能 。 该 研究 可 为 玉米 早期 干旱 、 营 养 胁迫 预警 提 
供 方 法 和 技术 支持 。 但 本 研究 不 足 在 于 光谱 数据 在 
室内 采集 ， 考 虑 到 田间 环境 往往 较为 复杂 ， 后 续 将 
对 模型 进行 迁移 学 习 以 提 高 其 在 不 同 环境 下 的 泛 化 
能 力 ， 力 争 开发 出 适用 于 智慧 农业 的 田间 干旱 、 营 
养 胁迫 在 线 监测 系统 。 


利益 冲突 声明 : 本 研究 不 存在 研究 者 以 及 与 公 
ae. 


参考 文献 : 
[1] BAR, KFA, NER, 等 . 近 红 外 高 光谱 成 像 扫 描 速 
度 对 拟 南 芥 TRA EKR I]. 光谱 学 与 光谱 分 
的 2020, 40(11): 3508-3514. 
LYU M Q, SONG Y J, WENG H Y, et al. Effect of near in- 
frared hyperspectral imaging scanning speed on prediction 


of water content in Arabidopsis[J]. Spectroscopy and spec- 

tral analysis, 2020, 40(11): 3508-3514. 
[2] 张 金 富 , Wowk, 王建 旭 , 等 . 基于 近 红 外 光谱 的 全 果 可 溶 
性 固形 物 特 征 波长 筛选 建 模 对 比分 析 [JMOL]. 激 光 与 兴 
电子 学 进展 . (2022-11-01)[2023-08-12]. https://kns.cnki. 
net/kcms/detail/3 1.1690.TN.20221031.1514.004.html 
ZHANG J F, TANG B, WANG J X, et al. Comparative 
analysis of characteristic wavelength screening of apple 


soluble solids based on near-infrared spectroscopy[J/OL]. 
Laser & Optoelectronics Progress. (2022-11-01) [2023-08- 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


THRI]. 作物 研究 , 2018, 32(4): 345-348. 
YANG J J, LU J W. Research progress of spectral technol- 
ogy in peanut quality detection[J]. Crop research, 2018, 32 


(4): 345-348. 
IMR, 金 夏 明 , EFF, 等 . 基于 高 光谱 图 像 的 生菜 叶片 
氮 素 含量 预测 模型 研究 [J]. 分 析 化 学 , 2014, 42(5): 
672-677. 


SUN J, JIN X M, MAO HP, et al. A model for predicting 
nitrogen of lettuce leaves based on hyperspectral imag- 
ing[J]. Chinese journal of analytical chemistry, 2014, 42 
(5): 672-677. 

GAO D H, LI M Z, ZHANG J Y, et al. Improvement of 
chlorophyll content estimation on maize leaf by vein re- 
moval in hyperspectral image[J]. Computers and electron- 
ics in agriculture, 2021, 184: ID 106077. 

TAN K Z, WANG S W, SONG Y Z, et al. Estimating ni- 
trogen status of rice canopy using hyperspectral reflec- 
tance combined with BPSO-SVR in cold region[J]. Che- 
mometrics and intelligent laboratory systems, 2018, 172: 
68-79. 

SHU M Y, DONG Q Z, FEI S P, et al. Improved estima- 
tion of canopy water status in maize using UAV-based dig- 
ital and hyperspectral images[J]. Computers and electron- 
ics in agriculture, 2022, 197: ID 106982. 
LICHTENTHALER H K, WELLBURN A R. Determina- 
tions of total carotenoids and chlorophylls a and b of leaf 
extracts in different solvents[J]. Biochemical society trans- 
actions, 1983, 11(5): 591-592. 

MASSAOUDI M, REFAAT S S, ABU-RUB H, et al. PLS- 
CNN-BiLSTM: An end-to-end algorithm-based savitzky- 
golay smoothing and evolution strategy for load forecast- 
ing[J]. Energies, 2020, 13(20): ID 5464. 

BAO Y D, KONG W W, LIU F, et al. Detection of glutam- 


ic acid in oilseed rape leaves using near infrared spectros- 


152 


智慧 农业 (中 英文 ) Smart Agriculture 


Vol. 5, No. 3 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


copy and the least squares-support vector machine[J]. In- 
ternational journal of molecular sciences, 2012, 13(11): 
14106-14114. 

CHEN Q S, JIANG P, ZHAO J W. Measurement of total 
flavone content in snow lotus (Saussurea involucrate) us- 
ing near infrared spectroscopy combined with interval 
PLS and genetic algorithm[J]. Spectrochimica acta part A: 
Molecular and biomolecular spectroscopy, 2010, 76(1): 
50-55. 

ARAÚJO M C U, SALDANHA T C B, GALVÃO RK H, 
et al. The successive projections algorithm for variable se- 
lection in spectroscopic multicomponent analysis[J]. Che- 
mometrics and intelligent laboratory systems, 2001, 57(2): 
65-73. 

ZARCO-TEJADA P J, BERJÓN A, LÓPEZ-LOZANO R, 
et al. Assessing vineyard condition with hyperspectral in- 
dices: Leaf and canopy reflectance simulation in a row- 
structured discontinuous canopy[J]. Remote sensing of en- 
vironment, 2005, 99(3): 271-287. 

MEHDAOUI R, ANANE M. Exploitation of the red-edge 
bands of sentinel 2 to improve the estimation of durum 
wheat yield in Grombalia region (Northeastern Tunisia) [J]. 
International journal of remote sensing, 2020, 41(23): 
8986-9008. 

GITELSON A A, KEYDAN G P, MERZLYAK M N. 
Three-band model for noninvasive estimation of chloro- 
phyll, carotenoids, and anthocyanin contents in higher 
plant leaves[J]. Geophysical research letters, 2006, 33 
(11): 431-433. 

FITZGERALD G, RODRIGUEZ D, O'LEARY G. Mea- 
suring and predicting canopy nitrogen nutrition in wheat 
using a spectral index: The canopy chlorophyll content in- 
dex (CCCI) [J]. Field crops research, 2010, 116(3): 
318-324. 

JORDAN C F. Derivation of leaf-area index from quality 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


of light on the forest floor[J]. Ecology, 1969, 50(4): 
663-666. 

DASH J, CURRAN P J. The MERIS terrestrial chloro- 
phyll index[J]. International journal of remote sensing, 
2004, 25(23): 5403-5413. 

GAO B C. NDWI: A normalized difference water index for 
remote sensing of vegetation liquid water from space[J]. 
Remote sensing of environment, 1996, 58(3): 257-266. 
HUNT E R, ROCK B N. Detection of changes in leaf wa- 
ter content using near- and middle-Infrared reflectances[J]. 
Remote sensing of environment, 1989, 30(1): 43-54. 
ZHANG J J, ZHANG W, XIONG S P, et al. Comparison 
of new hyperspectral index and machine learning models 
for prediction of winter wheat leaf water content[J]. Plant 
methods, 2021, 17(1): 1-14. 

PENUELAS J, PINOL J, OGAYA R, et al. Estimation of 
plant water concentration by the reflectance Water Index 
WI (R900/R970) [J]. International journal of remote sens- 
ing, 1997, 18(13): 2869-2875. 

ZARCO-TEJADA P J, USTIN S L. Modeling canopy wa- 
ter content for carbon estimates from MODIS data at land 
EOS validation sites[C]// IGARSS 2001. Scanning the 
Present and Resolving the Future. Proceedings. IEEE 
2001 International Geoscience and Remote Sensing Sym- 
posium (Cat. No. 01CH37217). Piscataway, New Jersey, 
USA: IEEE, 2002: 342-344. 

SUN H, FENG M C, XIAO L J, et al. Assessment of plant 
water status in winter wheat (Triticum aestivum L.) based 
on canopy spectral indices[J]. PLoS one, 2019, 14(6): ID 
e0216890. 

MENZE B.H., KELM B.M., MASUCH R. et al. A com- 
parison of random forest and its Gini importance with 
standard chemometric methods for the feature selection 
and classification of spectral data[J]. BMC Bioinformat- 
ics, 2009, 10(1): 1-16. 


A Hyperspectral Image-Based Method for Estimating Water and 
Chlorophyll Contents in Maize Leaves under Drought Stress 


WANG Jingyong', ZHANG Mingzhen , LING Huarong’, WANG Ziting”, GAI Jingyao™ 


(1. School of Mechanical Engineering, Guangxi University, Nanning 530004, China; 2. College of Agriculture, Guangxi Univer- 


Abstract: 


sity, Nanning 530004, China; 3. Guangxi Key Laboratory of Sugarcane Biology, Guangxi University, Nanning 530004, China ) 


[Objectives] Chlorophyll content and water content are key physiological indicators of crop growth, and their non-destructive detec- 


tion is a key technology to realize the monitoring of crop growth status such as drought stress. This study took maize as an object to 


develop a hyperspectral-based approach for the rapid and non-destructive acquisition of the leaf chlorophyll content and water content 
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for drought stress assessment. 

[Methods] Drought treatment experiments were carried out in a greenhouse of the College of Agriculture, Guangxi University. Maize 
plants were subjected to drought stress treatment at the seedling stage (four leaves). Four drought treatments were set up for normal 
water treatment [CK], mild drought [W1], moderate drought [W2], and severe drought [W3], respectively. Leaf samples were collect- 
ed at the 3", 6", and 9" days after drought treatments, and 288 leaf samples were collected in total, with the corresponding chlorophyll 
content and water content measured in a standard laboratory protocol. A pair of push-broom hyperspectral cameras were used to col- 
lect images of the 288 seedling maize leaf samples, and image processing techniques were used to extract the mean spectra of the leaf 
lamina part. The algorithm flow framework of "pre-processing - feature extraction - machine learning inversion" was adopted for pro- 
cessing the extracted spectral data. The effects of different pre-processing methods, feature wavelength extraction methods and ma- 
chine learning regression models were analyzed systematically on the prediction performance of chlorophyll content and water con- 
tent, respectively. Accordingly, the optimal chlorophyll content and water content inversion models were constructed. Firstly, 70% of 
the spectral data was randomly sampled and used as the training dataset for training the inversion model, whereas the remaining 30% 
was used as the testing dataset to evaluate the performance of the inversion model. Subsequently, the effects of different spectral pre- 
processing methods on the prediction performance of chlorophyll content and water content were compared. Different feature wave- 
lengths were extracted from the optimal pre-processed spectra using different algorithms, then their capabilities in preserve the infor- 
mation useful for the inversion of leaf chlorophyll content and water content were compared. Finally, the performances of different 
machine learning regression model were compared, and the optimal inversion model was constructed and used to visualize the chloro- 
phyll content and water content. Additionally, the construction of vegetation coefficients were explored for the inversion of chloro- 
phyll content and water content and evaluated their inversion ability. The performance evaluation indexes used include determination 
coefficient and root mean squared error (RMSE). 

[Results and Discussions] With the aggravation of stress, the reflectivity of leaves in the wavelength range of 400~1700 nm gradual- 
ly increased with the degree of drought stress. For the inversion of leaf chlorophyll content and water content, combining stepwise re- 
gression (SR) feature extraction with Stacking regression could obtain an optimal performance for chlorophyll content prediction, with 
an R° of 0.878 and an RMSE of 0.317 mg/g. Compared with the full-band stacking model, SR-Stacking not only improved R° by 2.9%, 
reduced RMSE by 0.0356mg/g, but also reduced the number of model input variables from 1301 to 9. Combining the successive pro- 
jection algorithm (SPA) feature extraction with Stacking regression could obtain the optimal performance for water content prediction, 
with an R° of 0.859 and RMSE of 3.75%. Compared with the full-band stacking model, SPA-Stacking not only increased R° by 0.2%, 
reduced RMSE by 0.03%, but also reduced the number of model input variables from 1301 to 16. As the newly constructed vegetation 
coefficients, normalized difference vegetation index(NDVI) [(Ry,)—Rss0)/(RyiotRss0)] and ratio index (RD (Ryoy/R,,7,) had the highest ac- 
curacy and were significantly higher than the traditional vegetation coefficients for chlorophyll content and water content inversion, re- 
spectively. Their R were 0.803 and 0.827, and their RMSE were 0.403 mg/g and 3.28%, respectively. The chlorophyll content and wa- 
ter content of leaves were visualized. The results showed that the physiological parameters of leaves could be visualized and the differ- 
ences of physiological parameters in different regions of the same leaves can be found more intuitively and in detail. 

[Conclusions] The inversion models and vegetation indices constructed based on hyperspectral information can achieve accurate and 
non-destructive measurement of chlorophyll content and water content in maize leaves. This study can provide a theoretical basis and 
technical support for real-time monitoring of corn growth status. Through the leaf spectral information, according to the optimal mod- 
el, the water content and chlorophyll content of each pixel of the hyperspectral image can be predicted, and the distribution of water 
content and chlorophyll content can be intuitively displayed by color. Because the field environment is more complex, transfer learn- 
ing will be carried out in future work to improve its generalization ability in different environments subsequently and strive to develop 


an online monitoring system for field drought and nutrient stress. 
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